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Abstract 

Zirconia (Zr02) is an important material with technological applications which are affected by point 
defect physics. Ab-initio calculations are performed to understand the structural and electronic properties 
of oxygen vacancies and interstitials in different charge states in cubic zirconia. We find oxygen interstitials 
in cubic Zr02 can have five different configurations - (110) dumbbell, (100) dumbbell, (100) crowd-ion, 
octahedral, and (111) distorted dumbbell. For a neutral and singly charged oxygen interstitial, the lowest 
energy configuration is the (110) dumbbell, while for a doubly charged oxygen interstitial the octahedral 
site is energetically the most favorable. Both the oxygen interstitial and the oxygen vacancy are negative-C/, 
so that the singly charged defects are unstable at any Fermi level. The thermodynamic stability of these 
defects are studied in terms of Fermi level, oxygen partial pressure and temperature. A method to determine 
the chemical potential of the system as a function of temperature and pressure is proposed. 
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1. Introduction 



Zirconia (Zr0 2 ) exists in three closely related phases with monoclinic, tetragonal and cubic 
symmetry. The monoclinic phase is thermodynamically stable near room temperature and up to 
1400 K, while at higher temperatures tetragonal and then cubic phases become stable. 1 Cubic 
zirconia is stable above 2570 K. 1 However, the cubic phase can be stabilized at much lower tem- 
peratures by doping Zr0 2 with divalent or trivalent cations like Y 3+ or Ca 2+ . 2 This cubic stabilized 
zirconia (CBZ) has high thermal shock resistance, high strength and toughness and has been un- 
der extensive investigation in recent years. 23 CBZ is also important for oxygen sensors and fuel 
cells. 4-6 To satisfy charge neutrality condition after cation doping, these materials contain oxygen 
vacancies and interstitials whose concentrations are related to the dopant concentration. 

Unstabilized cubic zirconia is also an interesting material for high temperature applications. A 
proposed methodology for exterior coatings in high speed aircraft involves coating with ZrB 2 /SiC 
ceramic. 7,8 Upon exposure to high temperatures in an oxidizing atmosphere, ZrB 2 /SiC oxidizes to 
Zr0 2 and Si0 2 +B 2 3 , a vitreous borosilicate glass which forms a sticky layer on the top. 9 Little 
is known about point defects and their stability for the high temperature phase of Zr0 2 , yet defect 
diffusion within the oxide scale is potentially a key issue for optimization of this material system, 
and experimental studies are difficult because of the extreme operating conditions. 

Moreover, the semiconductor industry has recently begun using hafnia and zirconia as a sub- 
stitute for silicon dioxide in MOSFET gate dielectrics, an application for which they have merit 
because of their high dielectric constant. 10-12 For this application the crystalline structure is gen- 
erally the low-temperature monoclinic phase. Also CBZ has good resistance to radiation induced 
damage, leading to a number of proposed and actual nuclear applications.— - — For both of these 
applications point defect physics are crucial. In the gate oxide application point defects can trap 
and scatter charge carriers and create other unfavorable effects. In the nuclear applications the 
rapid recombination of radiation-induced defects gives zirconia its radiation resistance. 

In this paper we study point defect formation in cubic Zr0 2 by using first principles density 
functional theory calculations. First we obtain the atomic structure and formation energy of the 
different intrinsic defects. Then we present results in which we calculate defect concentrations as a 
function of Fermi level, oxygen partial pressure and temperature, by studying the thermodynamic 
equilibria for the various defects with each other, with ambient oxygen, and with charge carriers. 
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2. Method of calculation 



We perform calculations using the Vienna Ab-initio Simulation Package (VASP),— r — a plane 
wave basis density functional theory (DFT) code, with generalized gradient approximation (GGA) 
and the exchange correlation potential of Perdew, Burke and Wang (PBE). 20 Spin polarization is 
taken into account in these calculations. The wave functions are expanded to a plane wave cut-off 
of 520 eV. Static relaxations are performed till the forces on individual atoms are smaller than 
10" 4 eV/A and total energy convergence within 10 5 eV is achieved. Initial k-point sampling 
is performed with 1 to 20 irreducible kpoints in the first Brillouin zone. The equilibrium lattice 
parameter of cubic Zr0 2 is 5.1276 A (with 20 irreducible kpoints) which is close to previous 
density functional calculations and experimental values.-^! All subsequent defect calculations are 
performed at this lattice constant, with no relaxation of the super cell unless otherwise mentioned. 
Total energy of a cubic Zr0 2 primitive cell converged to within 6 meV with 4 irreducible kpoints. 
Hence, we use 4 irreducible points for all subsequent calculations. 

The band gap is given by the energy difference between the highest occupied electronic energy 
level and the lowest unoccupied energy level. We obtain a band gap of 3.26 eV which is smaller 
than 5.40 eV obtained using GW calculations . 10112 ' 22 DFT calculations are known to underestimate 
the band gap. 12 To describe the DOS of the defect-free cell obtained from VASP calculations, a 
rigid shift of the conduction band edge is performed to obtain the experimentally reported band 
gap. The defect levels are not scaled while shifting the conduction band edge.^i 2 . 

Understanding the thermodynamic stability of charged defects involves incorporation or re- 
moval of electrons from the defect-free materials surrounding the defect. In the DFT calculations, 
the electronic charge neutrality of the supercell is achieved by a uniform jellium of equal and op- 
posite charge. This results in long range interaction between the defect and its periodic images 
in the calculation which must be subtracted off from the total energy of the supercell E vasp . Thus 
E = E vasp -E COI is the corrected energy of the system. -E vasp is the energy obtained from VASP 
calculations for a sample containing defects. For a single defect the correction term is 23 

q 2 a 

Aior — — "j j U) 

where q is the charge on the defect, a = 5.04 is the Madelung constant, L is the length of the 
cubic supercell used for calculations, k is the dielectric constant of cubic Zr0 2 and e Q is the per- 
mittivity of free space. For cubic zirconia, we take k = 36.8.- Typical Madelung corrections in 
our supercells are of the order of 0.05g 2 eV, where q is the charge of the defect (...,-2,-1,0,1,2,...). 



Given our primary interest in the high temperature cubic phase we calculate formation free 
energies for the lattice defects in zirconia as a function of temperature and oxygen partial pressure. 



3. Defect geometry, electronic charge state and stability 

Stability of point defects have to be analyzed with respect to some reference state. Since cubic 
zirconia is a high temperature phase and frequently used in contact with atmosphere we use oxygen 
gas molecule as the reference state for calculations of formation energies of oxygen vacancies and 
oxygen interstitials. 



3A. Reference state calculations 



First-principles calculations of an oxygen molecule are performed at T = K. The triplet state 
is found to be the ground state of O2. We find i?o a 2 sp =-9.86 eV and a bond length of 1.23 A. The 
singlet state is 1 .06 eV higher than the triplet state. 

The chemical potential (fio 2 ) of an oxygen gas molecule changes with temperature as well as 
the oxygen partial pressure. An estimate of the free energy of an oxygen molecule can be made 
based upon the quantum formulation of an ideal gas. The free energy of an ideal diatomic gas 
has contributions from the translational free energy, vibrational free energy and rotational free 
energy. 24 So the chemical potential of an oxygen molecule as a function of temperature (T) and 
pressure (P) is given by : 24 

hjJ 
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where, A is the thermal de-Broglie wavelength, k B is the Boltzmann constant, h is Planck's con- 
stant and I is the moment of inertia of the oxygen molecule. The chemical potential can also be 
expressed in terms of a reference pressure P r : 

P 

H = n r + k B T log — (3) 

where, fi r is the chemical potential at reference pressure. The reference pressure is generally taken 
as 1 atm. From ab-initio calculations, we find the vibrational frequency of 2 is u Q = 1568.06 
cm" 1 , which is close to experimental result of co Q =1565.4 cm" 1 . 25 For our subsequent analysis, 
the chemical potential of an oxygen atom is taken as half of fxo 2 : 

fio(T,P) = ^oAT,P) (4) 



Fig. [T] shows the variation in the oxygen chemical potential with temperature at different oxygen 
partial pressures. From this analysis we obtain fio(T = 300 K, P = latm) = -5.16 eV and no(T = 
2600 K, P = latm) = -8.20 eV. Since, the cubic phase in undoped Zr0 2 is stable above 2570 K, 
we have used T = 2600 K, for subsequent analysis. 

3B. Oxygen Vacancy 

Oxygen vacancies in cubic Zr0 2 have been discussed previously because their motion is related 
to the ionic conductivity of zirconia. They have also been studied to understand the defect levels, 
which are important for applications in the semiconductor industry In the low temperature 
monoclinic phase there are two non-equivalent oxygen lattice sites, and hence oxygen vacancies 
exhibit different coordination numbers.— In contrast, all oxygen vacancy sites in cubic zirconia 
are equivalent, and have four-fold coordination. 

In accord with previous work, we find that point defects significantly distort the nearby 
latticed— For different charge states the lattice relaxations around the vacancy are qualitatively 
different. Formation of neutral vacancy, V Q (using Kroger- Vink notation 29 ) results in inward dis- 
placement of the nearest neighbors. Allowing the volume of supercell to relax results in decrease 
in total volume by 1.03 A 3 . The neighboring Zr and oxygen ions move towards the oxygen va- 
cancy. In contrast to V D , the atomic structures of Yl + and V 2+ have outward movement of Zr ions 
while the oxygen ions move inwards. With increase in charge on the vacancy we find increase in 
distortion of the lattice. 

The relaxation volume of a charged vacancy has contributions from the local atomic relaxation 
around the vacant lattice site as well as the delocalized electrons present in the system. For a 
doubly charged oxygen vacancy the net change in volume is AVol(V 2+ ) = AVol(Zr0 2 +V 2+ ) + 
2AVol(Zr0 2 +e~). Here, AVol(Zr0 2 +e _ ) is the change in volume of the Zr0 2 cell with an ex- 
tra electron and AVol(Zr0 2 +V 2+ ) is the change in volume after incorporating a doubly charged 
vacancy (cell has an effective charge of +2e). This calculation yields AVol(V 2 1 + ) = 5.51 A 3 and 
AVol(V* + ) = 2.17A . The decrease in energy during this relaxation process is 0.09 eV for V„ + , 
0.32 eV for V* 4 " and 0.41 eV for V c (the cost of simultaneously putting electrons in the conduc- 
tion band is not included). The simulation cell following volume relaxation preserves the cubic 
symmetry. Thus vacancies tend to stabilize the cubic phase in Zr0 2 . 
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The process of vacancy formation (say Vl + ) can be described by the following equation : 

o ^i0 2 ( ga s)+V2+ + 2e- (5) 

Physically, this means the atmosphere acts as sink for the oxygen atoms removed from Zr0 2 - The 
removal of a neutral oxygen atom from the system (leaving behind a neutral Zr0 2 lattice) results 
in a vacant lattice site. The vacant lattice site can trap electrons to form a neutral vacancy or singly 
charged vacancy. 

Fig. |2(a)| and Fig. |2(b)| shows the formation energy of oxygen vacancy as a function of Fermi 
level for the charge states V G , V* + , at T = K and T = 2600 K. When the Fermi level is 
near the valence band edge, is the most stable charge state of a vacancy, but at higher Fermi 
levels neutral charge state (V ) becomes energetically favored. We define e (V^/V^) as the 
point at which and V 1+ have similar values of formation energy. It signifies the transition in 
thermodynamic stability of doubly charged vacancy to a singly charged vacancy with increase in 
Fermi energy. Also from Fig. |2(a)| and Fig. |2(b)| we can see that the thermodynamic transition 
level 30 e (V„ + /Vo + ) is located higher than e (V^/V ). Hence, as in the monoclinic phase, the 
formation of is not favorable in cubic zirconia at any Fermi level, meaning that the oxygen 
vacancy is a negative-f/ defect. 10 The defect levels we find are similar to previous estimates using 
screened exchange method and weighted density approximation.— 



3C. Oxygen Interstitial 

Cubic Zr0 2 has CaF 2 structure with metal atoms in a face centered cubic sub-lattice and oxygen 
atoms occupying a simple cubic sub-lattice. Since the oxygen sub-lattice is not close packed, 
oxygen defects, rather than zirconium defects, play a dominant role in determining the properties 
of this oxide. We find that oxygen interstitials in cubic zirconia have different possible atomic 
structures than those previous reported for the monoclinic phase. We have studied different oxygen 



interstitial configurations and geometries : (i) (110) dumbbell (Fig. |3(a)[ ), (ii) (100) dumbbell 



(Fig. |3(b)| ), (Hi) (100) crowd-ion (Fig. |3(c)[ ) and (iv) octahedral site (Fig. |3(d)[ ). For ease of 



reference the interstitial configurations are labeled as O/n ), O(i o) respectively for the (110), and 
(100) dumbbell configurations, O(i o c ) for (100) crowd- ion and O oct for octahedral site. Another 
possible interstitial configuration is (111) dumbbell. While (110) dumbbell and (100) dumbbell 
have dumbbell centers on an original atomic site, a (111) dumbbell configuration upon relaxation 
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has dumbbell center shifted and not on an original atomic site. (Ill) dumbbell in fact looks quite 
similar to oct , but not equivalent. It has 0.88 eV higher energy than O/no) and the dumbbell 
separation is 1.43 A 3 . 

The relative stability of different interstitial configurations with respect to Omo) is compared 
in TablelH For an oxygen interstitial with neutral charge, Oni ) configuration is energetically most 
favorable. The dumbbell separation is close to the 0-0 bond length in peroxides Incorporation 
of an extra oxygen atom results in outward movement of the nearest neighbor Zr ions. In this case 
the lattice distortion is symmetrical about (110) plane. The Zr ions along (110) moves 6% (of the 
Zr-0 bond length in an ideal crystal) outward. The nearest neighbor oxygen ions move slightly 
inward by 0.7%. The two oxygen ions forming the dumbbell structure are separated by 1.46 A. 
The change in volume of supercell for Oni ) was calculated by allowing the supercell to relax. 
For O(no) structure, we find it to be 36.10 A 3 , much higher than the values obtained for vacancies. 
Volume relaxation destroys the cubic symmetry of the simulation cell. 

The electronic local density of states plots near one of the dumbbell ions are the same for spin 
up and spin down. This is also true for any oxygen atom in the bulk far away from the interstitial 
configuration. Thus the oxygen ions in the system do not have residual spins. Robertson et al. 
performed calculations for oxygen interstitials in cubic zirconia and reported electronic properties 
of O(ioo), though our findings suggest the O(no) as the most stable. 

In order to verify the lowest energy configuration we perform ab-initio molecular dynamics 
(MD) simulations at T = 2000 K for 1.5 ps. The high temperature enhances mobility of the 
ions and the structures obtained are much different from those at T = K. Next, we select a few 
intermediate configurations from the MD simulation results and cool them to T = K in molecular 
dynamics and perform static energy minimization. This procedure typically yields significantly 
lower-energy structures (AE ~ 0.5 eV) for vacancies and interstitials. These low energy structures 
have low symmetry (due to Jahn-Teller distortion). This methodology is well-behaved in that when 
the same procedure is applied to a perfect cubic zirconia cell, the perfect cell is recovered. When 
performing this procedure on defects we do not allow the cell lattice vectors to relax, since if this 
is allowed the cell can become tetragonal. 

For singly charged interstitials, we find that the O| no ^ configuration is energetically most fa- 
vorable in contrast to O^ reported previously by Robertson et al. 12 We come to this conclusion 
by performing MD simulations at T = 2000 K and then quenching to T = K. After static relax- 
ation, we find the separation between the ions in a dumbbell is 1 .98 A, much higher than that in a 
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Defect 


Interstitial 


0<iio) 


0(m> 


O(ioo) 


O(100c) 


Ooct 


o? 


Energy (relative, eV) 





0.88 


0.95 


1.92 


3.02 




Dumbbell separation (A) 


1.46 


1.44 


1.43 








Energy (relative, eV) 







3.80 


3.79 


1.90 




Dumbbell separation (A) 


1.98 




1.93 







TABLE I: Relative stability and dumbbell separations for neutral and singly charged oxygen interstitials 
in cubic zirconia. After annealing by ab-initio MD and possessing lower symmetry due to Jahn-Teller 
distortion. 

neutral cell. The 0* m ^ configuration becomes unstable and the interstitial atom moves to nearest 
octahedral site. Table U shows the relative stability of different configurations for a singly charged 
interstitial with respect to O* no ^ . 

The most stable site for is the octahedral site. Interstitial atoms placed at the other config- 
urations become unstable and move to occupy the nearest octahedral site upon relaxation. Similar 
to the neutral oxygen interstitial, the projected density of states for spin-up and spin-down are 
identical and hence there is no residual spin on the interstitial atom. Similar to the case of charged 
vacancies, incorporating charged interstitials, Of~ results in large distortion of the parent lattice. 
The neighboring oxygen and zirconium ions display outward distortions. 

Physically, the process of formation of Of ~ can be understood by this reaction : 

io 2 (gas)^0^ + 2h 1+ . (6) 

Here, the atmosphere acts as source of oxygen. The system has to be neutral before and after the 
formation of the interstitial, so formation of a charged interstitial results in the formation of holes. 
Fig. |2(a)| and Fig. |2(b)| shows the formation energy of 0?~, O; 1 and O? as a function of the Fermi 
energy of the system at 1 atm pressure and T = K and T = 2600 K respectively. With increasing 
temperature incorporation of oxygen interstitials become energetically expensive. For example, 
the incorporation of neutral interstitial at T = K costs E = 0.86 eV but at T = 2600 K is E = 4. 12 
eV. From the formation energy plot, we conclude that the thermodynamic defect level e (O? /Ol~) 
at T = K is higher than e (O^/O? - ) . Hence, as in the monoclinic phased, 0* is energetically 
not favored at any Fermi level in cubic Zr0 2 , meaning that the oxygen interstitial is a negative-?/ 
defect. 

Table |n] shows some of the possible reactions between the anionic defects in Zr0 2 and their 

8 



No. 


Reaction 


Energy (eV) 


1 


0?-+OP=^20, 1- 

11 1 


0.83 


2 


V2++V^2Vi+ 

o o o 


0.43 


3 


Null =>0°+V° 

1 o 


E''=5A5 


4 


Null = > 0}-+V 1 + 

1 o 


££=4.08 


5 


Null =^0?~+V;? + 

1 o 


£7p=1.75 


6 


0°+V°^0 1 1 ~+Vi+ 


-1.07 


7 


Oj 1_ +Vj+ =>o 1 2 ~+v 2+ 


-2.33 


8 


o°+v° =^o 2 -+v 2 + 

w l O o 


-3.40 


9 


Null =^e~+h+ 


E g =5A 



TABLE II: Defect reactions and the corresponding energy change at T = K. The energy change is obtained 
from the formation energy of the defects. The values reported have been corrected for dipole interactions. 

reaction energies. These reaction energies are calculated at K and zero pressure assuming dilute 
concentration of the defect species. A positive energy means the forward reaction is unfavor- 
able. Reactions 1 and 2 have positive reaction energies, which means that the formation of singly 
charged interstitials and vacancies are not favored. Thus, two singly charged interstitials or va- 
cancies would decay into neutral and doubly charged defects. Reactions 6-8 suggest that a pair of 
neutral vacancy and interstitial would prefer to form a pair of doubly charged Frenkel defect. 

Next we study the incorporation of an oxygen molecule. The formation energy of an oxygen 
molecule in cubic Zr0 2 is 10.42 eV at 1 atm oxygen partial pressure and T = 2600 K. At zero 
pressure and T = K the formation energy is 4.5 1 eV, much higher than interstitial formation 
energy of 0.86 eV in a neutral cell. The oxygen molecule formation energy in cubic Zr0 2 is 
close to the values reported for low temperature monoclinic phase of Hf0 2 and much higher than 
those reported to silica. 33 The bond length of the oxygen molecule is 1.27 A. By performing MD 
simulations at 2600K for lps and then subsequent relaxation, we find that the molecule prefers to 
dissociate into two (110) dumbbells. The corresponding decrease in energy obtained from static 
relaxation is 1.93 eV. 
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4. Defect Equilibrium Reactions 



The thermodynamic stability of different point defects can be understood by writing down the 
equilibrium chemical reactions for the formation of these defects. Analysis of the defect chemistry 
is done by obtaining the laws of mass action for different competing processes.- 3 -^— Apart from 
the laws of mass action, the various point defects present must also satisfy the electronic charge 
neutrality condition. Here we consider the stable anionic defects in Zr0 2 , namely O i5 2 ~, V Q and 
V 2+ . Each defect is considered to be in its lowest energy configuration. In cubic Zr0 2 , the metal 
ions occupy the closed packed face centered lattice sites. So cationic defect formation energies 
are much higher than anionic defect formation energies. 37 Hence we neglect the cationic defect 
contributions in subsequent analysis. 

Possible reactions involving O 2- and V 2+ are as follows: 

• formation of Frankel defects in the anionic sub-lattice : 

o Of- + V 2+ K F = [ ° i2 "] v [ 2 V ° 2+] (7) 

• formation of an oxygen vacancy by losing oxygen to the atmosphere : 

<=^0 2 (ga8)+Vr+2e- K, = V , 2 W 

• formation of electron-hole pair : 



Null <{=► e" + h + Ki= [h J t ] = e~ E ^ T (9) 



N V N C 



• formation of oxygen interstitial : 

-0 2 (gas)^Or + 2h+ K2 = ^K^f 00) 

• and charge neutrality condition : 

2[0 ; ] 2 -+ [e-] = 2 [V ] 2+ + [h+] (11) 

where, [e~] and [h + ] are the electron and hole concentration, respectively, [Oi 2- ] is the oxygen 
interstitial concentration, [V D 2+ ] is the oxygen vacancy concentration, iV a is the concentration of 
anionic sites and E s is the band gap. N c = 2 (27rm* c k B T/h 2 f /2 and N v = 2 (2irmlk B T / h 2 f' 2 
are the constants used to normalize the electron and hole concentration. 36 
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Since K F = e ~ El? / kBT and K . _ e -E g /k B T^ usmg va i ues f rom Table HH we find K F > K,. 
Hence the concentration of electronic defects [h + ], [e _ ] are much smaller than the doubly charged 
Frankel defects. In this limit, at intermediate oxygen partial pressure, the charge neutrality equa- 
tion can be simplified (Brouwer approximation 38 ) to [0 ; 2 ~] = [V Q 2+ ] oc \/K F . When the oxygen 
partial pressure is very high, (flOl) shifts towards right and we get more oxygen interstitials. So us- 
ing the charge neutrality (JTTJ) condition, the concentration of oxygen interstitials can be expressed 
in terms of the oxygen partial pressure: 

[V 2 +] oc P 02 1/6 [C, 2 "] oc P$. (12) 

The defect concentrations are dependent on multiple factors. In nature, the defect concentra- 
tions can be changed by adding dopant or by varying the external environment (temperature and 
pressure). Here we take undoped zirconia and assume that the defect concentrations, are deter- 
mined by these two external degrees of freedom. 

For the neutral defects the thermodynamic equilibrium reactions for O; and V G are given by, 

• formation of Frankel defects in the anionic sub-lattice : 

o ^o i + v (13) 

• formation of an oxygen vacancy by loosing oxygen to the atmosphere : 

o <=► -0 2 (gas) + V K'; = 1 oi °> (14) 

• formation of oxygen interstitial : 

l -0 2 (gas) <=► O, Kl = J^L (15) 

2 NJ>£ 

Analyzing the system for neutral defects is much simpler than the charged defects. Solving the 
above we get 

[Oi] oc Ftf [V ] oc P l J 2 (16) 

5. Determination of equilibrium Fermi energy 

The formation energy plots provide information about the stability of different point defects as 
a function of temperature, pressure and Fermi energy of the system. Consider for example, defect 
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formation energies at 2600 K and 1 atm (Fig. |2(b)[ ). The doubly charged defects have the lowest 
energies. However, for Fermi energy values less than ~2.1 eV and more than 3 eV the system 
becomes unstable leading to spontaneous formation of defects. This pins the Fermi energy of the 
system to a small region of 2.1 to 3 eV. The overall system needs to maintain charge neutrality. 
For this to happen the Fermi energy (/x e ) of the system needs to be 2.6 eV so that doubly charged 
anionic defects have same concentration. 

In order to understand this better let us look at the charge neutrality condition for undoped Zr0 2 
described in (fTTI) . We neglect the singly charged defects as they are thermodynamically not stable. 
We do not include the cationic defects as higher formation energy would make their concentration 
insignificant. As discussed in before, we can have three regimes : (a) [e~] = 2 [V Q ] 2+ - this 
implies significant electronic conductivity which is not physical for an insulator like Zr0 2 ; (b) 
[h + ] = 2 [Oi] 2 - similar to the above condition this is not possible; (c) [O;] 2- = [V Q ] 2+ - this 
regime seems physically plausible. 

In order to confirm this conclusion, we plot the defect concentrations at 2600 K and 1 atm 
pressure as a function of the Fermi energy. Fig. [5] shows a comparison between the defect concen- 
trations and the electron and hole concentrations. The defect concentrations are calculated from 
the rate constants derived from the laws of mass action, namely, ® and (flOl) . The electron and 
hole concentrations are obtained by [e~] = N c e~( Ec ~^ /kBT and [h + ] = N v e~^~ Ec ^ kBT '. Clearly, 
the doubly charged defects have much higher concentration than the neutral and singly charged de- 
fects. As discussed above, this analysis shows that the electron and hole concentrations are much 
smaller than the anionic defect concentrations and hence conductivity is by migration of ions. 

6. Conclusions 

An ab-initio study of anionic point defects in cubic zirconia are reported. For O; and Ol~ 
the energetically favorable configuration is (110) dumbbell and for Of - it is an octahedral site. 
Using a combination of high temperature molecular dynamics and static minimization we obtain 
the global minima of the potential energy for these defects. Our ground state configurations differ 
substantially from those previously reported by Robertson et al. 12 

Oxygen interstitials and vacancies in the cubic phase exhibit negative-?/ behavior, meaning 
that the singly charged defects are suppressed at any Fermi level. From the analysis of possible 
reactions between these anionic defects we find that singly charged defect species are thermo- 
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dynamically less preferable than a combination of neutral and doubly charged defects. Using 
molecular O2 as the reference state, defect formation energies are obtained at K and 2600 K. 
The neutral interstitials have ~ 2 eV higher formation energy than the vacancies at 2600K and 1 
atm pressure. The system can lower its energy if the oxygen ions forming the interstitial dumbbell 
diffuse to the surface and move out to the atmosphere as neutral oxygen molecule. This results 
in the formation of a neutral vacancy. The corresponding change in energy is AE = E(Zr02+V o ) 
- [E(Zr02+Oj) - E(02)] = -2.49 eV. Hence, we conclude that neutral interstitials are not stable in 
cubic zirconia. 

We find the relaxation volumes for the anionic vacancies are much smaller than the oxygen 
interstitials. Unlike the oxygen interstitials, the local atomic relaxation around oxygen vacancies 
tend to stabilize the cubic phase. We find a charged vacancy has considerably higher relaxation 
volume compared to a neutral vacancy. Higher relaxation volume signifies that ZrC>2 is susceptible 
to diffusional creep, which can degrade its mechanical properties with increase in temperature. 

Zirconium based materials are attractive for nuclear waste disposal. The Zr-0 binary phase 
diagram 39 shows the existence of non- stoichiometric Zr0 2 -s phase suggesting vacancies are dom- 
inant defect species. Assuming vacancies are dominant defects, we find that change in volume 
(AV) due to presence of defects in cubic Zr0 2 is not large (AV/V<0.5%). Similar results have 
also been obtained for ZrSi04.— 

For a stable system, formation energy of a defect has to be positive. This puts a constrain on 
the accessible Fermi energy values. From Fig. [2] we see that at K, the feasible Fermi energy 
values are ~0.5-l .5 eV, while at 2600 K Fermi energy of the system lies in the range of ~2. 1-3.0 
eV. For undoped ZrC>2, the Fermi energy would lie mid-gap, however, the Fermi energy values at 
low temperature suggest that concentration of doubly charged defects can be tuned by doping with 
acceptor impurities. 

Equilibrium defect concentrations in undoped ZrC>2 are obtained as a function of oxygen partial 
pressure and temperature by solving the laws of mass action. Based on the charge neutrality 
condition a method to calculate the Fermi energy of the system is proposed. It is important to note 
that these defect concentrations are evaluated without accounting for defect interactions. Thus, in 
the high concentration limit, the actual concentrations can be different from those shown in the 
figure. 

Our results are also relevant to cubic stabilized zirconia (CBZ), in which a p-type dopant such 
as Y or Ca lower the Fermi level. These low values of Fermi level signifies the abundance of dou- 
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bly charged oxygen vacancies. Thus, our results demonstrate that, by solving defect equilibrium 
equation separately and then comparing their concentrations, an estimate of the dominant defect 
species can be made. 
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FIG. 1 : The chemical potential of an oxygen atom as a function of temperature at different oxygen partial 
pressures calculated from the quantum formulation of an ideal gas. (fio = 0.5//o 2 ) 
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(a) 



(b) 




FIG. 3: Snapshots of oxygen interstitial configurations in cubic zirconia viewed along (100) direction. 3(a) 



The (110) dumb bell configuration neutral oxygen interstitial. This is the lowest energy configuration of the 
neutral oxygen interstitial, which is energetically favored at low Fermi level. |3(b)| (100) crowdion configu- 
rations of neutral oxygen interstitial. 3(c)| (100) crowdion configurations of neutral oxygen interstitial. |3(d)| 
Doubly charged oxygen interstitial in octahedral configuration. This is the ground state of a doubly charged 
oxygen interstitial. 
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FIG. 4: The chemical potential (Fermi energy) of an undoped cubic Zr02 system as a function of 
ture at different oxygen partial pressures. 
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FIG. 5: Comparison of point defect concentration as a function of Fermi level of the system at 1 atm 
pressure and[5]T = 2600 K. The concentration profiles are obtained by solving each set of defect equilibrium 
equations separately. The doubly charged defects : (solid red line) and (dotted red line) have 
higher concentration than the singly charged (V^ + - solid blue line, and O; 1 - dotted blue line)and neutral 
defects (V Q ) - slid green line, and Oi - dotted green line. The defect concentrations are expressed as fraction 
of anionic sites in the system. 
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